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The  variation  in  the  average  particle  number  and  energy  flux  in  the  high-latitude  region  is  determined 
as  a  function  of  the  solar  wind  velocity,  Viw,  and  the  z  component  of  the  interplanetary  magnetic  field, 
B;.  The  study  is  made  using  the  data  from  the  SSJ/4  detectors  on  the  satellites  of  the  Defense 
Meteorological  Satellite  Program.  In  the  study  the  high-latitude  region  is  divided  into  spatial  elements 
in  magnetic  local  time  and  corrected  geomagnetic  latitude.  One  such  matrix  of  spatial  divisions  is 
constructed  for  each  of  30  paired  ranges  of  values  in  Vsw  andB..  There  are  six  divisions  in  B.  covering 
the  range  from  - 10  to  + 10  nT  and  five  divisions  in  velocity  covering  the  range  from  200  to  800  km/s. 
Using  approximately  34  million  SSJ/4  spectra,  the  average  electron  and  ion  spectra  are  determined  in 
each  spatial  element  for  each  of  the  30  paired  values  of  Vsw  and  B,.  From  the  average  spectra,  the 
average  integral  energy  and  number  fluxes  for  electrons  and  ions  are  calculated  in  each  spatial  bin. 
These  values  are  then  spatially  integrated  to  give  the  average  hemispheric  inputs  of  the  integral  particle 
energy  and  number  flux.  Both  quantities  are  found  to  vary  in  a  simple  and  consistent  manner  with  both 
Vsw  and  Bj.  For  both  electrons  and  ions  the  variation  with  B.  tends  to  a  minimum  value  for  weak  to 
moderately  strong  B.  positive.  For  decreasing  values  of  B.  from  the  minimum,  both  quantities 
increase  at  a  rate  greater  than  linear.  For  increasing  values  of  B.  from  the  minimum,  both  quantities 
tend  to  increase,  least  for  the  hemispheric  electron  energy  flux  and  most  for  the  hemispheric  ion 
number  flux.  The  variation  in  both  quantities  with  for  both  electrons  and  ions  is  generally  linear 
with  the  steepest  slopes  for  the  electron  hemispheric  energy  flux.  The  variation  with  B,  and  V5W  for 
all  four  quantities  can  be  well  fit  by  a  simple  quadratic  equation  either  of  the  form/fB;,  Vsw)  =  n(B. 
-  b)‘  +  cVsw  +  d  or  /IB,,  Vsw)  =  l«(B.  -  b )2  +  IHcVjw  +d].  Either  form  indicates  that  the 
variation  is  no  greater  than  quadratic  in  B.  and  no  more  than  linear  in  Vsw  and  that  the  variation  is 
approximately  symmetric  for  values  above  and  below  the  minimum  value  of  B . .  namely,  b .  This  result 
is  significantly  different  than  either  the  half-wave  rectifier  or  epsilon  function  for  describing  the 
solar-wind  magnetospheric  interaction. 


1.  Introduction 

This  study  is  designed  to  quantify  the  response  of  auroral 
precipitation  to  varying  solar  wind  conditions.  We  parame¬ 
terize  auroral  precipitation  by  the  precipitating  electron  and 
ion  hemispheric  energy  flux  and  number  flux  and  the  solar 
wind  condition  by  the  z  component  of  the  interplanetary 
magnetic  field  (IMF  Bz)  and  the  solar  wind  velocity  (Vsw). 
We  have  previously  quantified  the  response  of  auroral 
precipitation  to  varying  geomagnetic  activity  using  the  Kp 
index.  This  work  may  be  found  in  the  papers  by  Hardy  et  al. 
[1985]  (electron  response)  and  Hardy  et  al.  [1989]  (ion 
response),  which  are  referred  to  as  paper  1  (electrons)  and 
paper  2  (ions)  throughout  the  remainder  of  this  paper. 
Although  sorting  by  Kp  provides  insight  into  the  dependence 
of  auroral  fluxes  and  morphology  on  geomagnetic  activity,  it 
does  not  directly  address  the  question  of  which  external 
mechanisms  are  responsible  for  driving  the  magnetospheric 
processes.  We  therefore  extend  our  work,  replacing  Kp  with 
B.  and  Fsw  as  ordering  parameters. 

It  is  widely  acknowledged  that  the  polarity  and  magnitude 
of  IMF  B.  play  the  major  role  in  controlling  the  merging 
process  between  closed  geomagnetic  held  lines  and  open 
interplanetary  magnetic  field  lines  [ Russell ,  1986;  Battmjo- 
hann  and  Haerendel.  1985].  As  such,  it  significantly  affects 
the  Earth’s  magnetic  field  topology  and  plays  a  fundamental 
role  in  modulating  the  strength  of  the  solar  wind- 
magnetosphere  coupling.  The  solar  wind  speed  also  plays  a 
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significant  role  both  electrodynamically,  through  the  im¬ 
pressed  electric  field  (VSWB:)  [Maezawa  and  Mttrayama, 
1986],  and  magnetohydrodynamically  [Spreiter  and  Alksne, 
1969]. 

Gussenhoven  [1988]  defined  the  ground  state  of  the  mag¬ 
netosphere  as  the  condition  resulting  from  a  period  of 
extended  minimal  solar  wind-magnetosphere  coupling  dur¬ 
ing  which  its  energy  content  has  decayed  to  a  minimum.  By 
examining  low-altitude  signatures  of  magnetospheric  activity 
(relating  to  both  auroral  oval  and  polar  cap  phenomena) 
under  variable  solar  wind  conditions,  she  concludes  that  the 
baseline  magnetosphere  occurs  after  a  period  of  low  B, 
north  (near  B.  =  0)  and  low  V,w.  The  polarity  of  B.  is 
critical  in  determining  the  manner  in  which  the  magneto¬ 
sphere  is  energized.  A  southward  B.  leads  to  an  energetic 
oval  and  a  quiet  polar  cap,  whereas  a  northward  B.  leads  to 
a  quiet  oval  and  a  more  active  polar  region.  Low-altitude 
signatures  of  magnetospheric  processes  driven  by  B.  south 
(active  auroral  oval)  and  by  strong  B.  north  (active  polar 
cap)  are  beginning  to  define  the  full  range  of  processes  by 
which  the  magnetosphere  is  energized.  An  accurate  assess¬ 
ment  of  the  dynamical  state  of  the  magnetosphere  must 
involve  consideration  of  activity  in  the  polar  cap  as  well  as  in 
the  auroral  oval. 

The  solar  wind  and  IMF  energize  the  magnetosphere 
through  a  number  of  interrelated  processes  which  include 
impressing  electric  potentials,  driving  currents,  and  energiz¬ 
ing  particles.  Despite  the  complexity  of  the  coupling,  high 
correlations  have  been  obtained  between  ionospheric  signa¬ 
tures  of  these  processes  and  solar  wind  parameters.  Low- 
altitude  measures  of  magnetospheric  activity  that  have  been 
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correlated  include  (1)  the  cross-polar  cap  potential,  (2)  the 
auroral  electrojet  current  indices  AE  (AL,  AU),  (3)  the 
field-aligned  current  intensity,  (4)  the  frequency  and  inten¬ 
sity  of  polar  cap  arcs,  and  (5)  the  precipitating  particles 
(auroral  boundaries  and  hemispheric  fluxes). 

Reiff  and  Luhmann  [1986J  have  reviewed  efforts  to  model 
the  cross-polar  cap  voltage  (Fpc)  as  a  function  of  B .  and 
Vsw.  They  report  that  for  Bt  south,  dayside  merging  is 
responsible  for  between  80  and  95%  of  Fpc  with  the  remain¬ 
der  due  to  nightside  magnetotail-lobe  merging  and  magneto¬ 
pause  boundary  layer  quasi-viscous  interaction.  For  B, 
north,  the  contribution  to  Vpc  from  dayside  merging  is 
greatly  diminished,  while  that  from  the  latter  two  processes 
increases  slightly,  so  that  functions  which  best  model  Vp,. 
are  those  that  do  not  cut  off  sharply  for  weak  B .  north. 

The  AL  and  AE  indices  have  been  shown  to  be  strongly 
correlated  to  Bz  and  [Maezawa,  1979;  Murayama,  1982; 
Baker,  1986,  and  references  therein).  As  substorm  indices, 
AE  and  AL  are  most  valuable  during  periods  of  B.  south 
when  the  enhanced  auroral  electrojet  currents  are  in  close 
proximity  to  the  ground  stations  measuring  their  correspond¬ 
ing  magnetic  perturbations.  However,  during  periods  of  B , 
north,  the  electrojets  diminish  and  the  region  of  enhanced 
magnetic  activity  typically  moves  poleward  beyond  the 
effective  range  of  the  ground  stations.  Thus  if  one  relies 
solely  on  the  AE  index  as  the  monitor  of  solar  wind- 
magnetospheric  coupling,  one  is  lead  to  model  the  magneto¬ 
sphere  as  a  half-wave  rectifier  where  energy  coupling  pro¬ 
ceeds  for  B,  south  but  turns  off  for  Bz  north. 

Field-aligned  current  systems  also  indicate  a  strong  depen¬ 
dence  on  IMF  B..  It  has  been  shown  [Zanetti  and Potemra, 
1986)  that  the  integrated  current  intensity  of  the  dayside 
auroral  region  1  current  system  maximizes  for  strong  B, 
negative  and  decreases  to  a  minimum  for  weak  B.  positive, 
at  which  point  the  "NBZ”  field-aligned  currents  develop 
poleward  of  the  auroral  oval  and  increase  with  increasing 
positive  B,.  These  NBZ  currents  characterize  the  polar  cap 
arc  state,  and  they  can  be  extremely  large  at  times  [lijima 
and  Shibaji,  1987;  Rich  et  al.,  1990).  Gussenhoven  [1982]  has 
shown  that  polar  cap  arcs,  as  viewed  from  auroral  satellite 
images,  occur  preferentially  (and  with  greater  intensity) 
during  periods  of  Bz  north  and  large  V'sw. 

The  variation  of  auroral  boundaries  with  solar  wind  pa¬ 
rameters  also  provides  information  on  the  effects  of  the  solar 
wind  on  magnetosphere  dynamics.  Precipitating  electrons 
have  been  used  to  define  various  auroral  boundaries  [Hardy 
et  al.,  1982]  which  map  to  dynamically  distinct  magneto- 
spheric  regions:  the  equatorward  boundary  maps  to  the  inner 
edge  of  the  central  plasma  sheet,  the  transition  boundary 
maps  to  the  transition  between  the  plasma  sheet  and  the 
boundary  layer,  and  the  poleward  boundary  maps  to  the 
furthest  extent  in  the  magnetotail  lobes  that  accelerated 
boundary  layer  electrons  are  found.  Hardy  et  al.  [1981] 
examined  the  relation  between  the  electron  equatorward 
boundary  and  the  solar  wind  parameters  (B.,  VWB.)  and 
found  an  expanding  oval  for  increasingly  negative  B:  but  no 
systematic  variation  for  Bz  >  I  nT.  Makita  et  al.  [1983] 
found  that  all  three  boundaries  systematically  moved  equa¬ 
torward  for  an  increasing  magnitude  of  Bt  negative.  How¬ 
ever,  for  an  increasing  magnitude  of  B.  positive,  only  the 
movement  of  the  poleward  boundary  toward  higher  latitudes 
was  moderately  well  correlated.  Makita  et  al.  [1988]  f  »jnd 
that  soft  electron  precipitation  (i.e.,  average  energy  of  less 


than  0.5-1  keV)  mapping  to  boundary  plasma  regions  is 
enhanced  for  periods  of  strong  Bz  north  relative  to  weak  Bz 
north. 

The  total  hemispheric  energy  and  number  flux  of  precipi¬ 
tating  electrons  and  ions  is  an  important  low-altitude  mea¬ 
sure  of  magnetospheric  activity  in  that  by  its  global  nature  it 
is  biased  toward  neither  the  active  oval  (B,  south)  nor  the 
active  polar  cap  ( B .  north).  Foster  et  al.  [1986]  have  statis¬ 
tically  determined  precipitating  electron  hemispheric  energy 
flux  as  a  function  of  Bz.  However,  they  did  not  consider  Fsw 
dependence  and,  owing  to  their  two  B,  bin  intervals  ( B ,  < 
-3  nT;  B,  >  +3  nT),  were  unable  to  test  for  minimum 
precipitation  levels  at  low  B t  values.  By  including  Vsw  as  a 
sorting  parameter,  increasing  the  B,  interval  to  include 
lower  absolute  magnitudes  of  B.  and  analyzing  ions  as  well 
as  electrons,  we  extend  the  work  of  Foster  and  coworkers  in 
using  auroral  precipitation  as  an  important  probe  of  mag¬ 
netospheric  dynamics. 

This  paper  is  divided  into  seven  sections.  Following  the 
introduction  (section  1),  we  present  ir  section  2  the  instru¬ 
ment  and  data  coverage;  in  section  3,  the  data  binning 
procedures;  and  in  section  4,  the  derivation  of  hemispheric 
fluxes  from  the  average  spectra.  The  results  are  presented  in 
section  5  and  discussed  in  section  6.  Section  7  is  the 
summary  and  conclusion. 

2.  Instrumentation  and  Data  Coverage 

The  data  analyzed  in  this  study  are  the  particle  fluxes 
measured  by  the  SSJ/4  ion  and  electron  detectors  flown  on 
the  F6  and  F7  satellites  of  the  Defense  Meteorological 
Satellite  Program  (DMSP).  The  time  period  of  F6  data 
coverage  includes  January  1983  to  July  1987,  and  for  F7, 
December  1983  to  July  1987.  The  DMSP  satellites  are 
continuously  maintained  in  circular,  polar,  sun-synchronous 
orbits  at  an  altitude  of  840  km  and  an  inclination  of  98.3°. 
One  orbit  (F6)  lies  in  the  0640-1840  magnetic  local  time 
(MLT)  meridian  plane,  and  the  other  (F7)  lies  in  the  1030— 
2230  MLT  meridian  plane.  The  satellites  are  three-axis 
stabilized  with  the  SSJ/4  sensor  apertures  oriented  toward 
the  local  zenith  at  all  times  to  measure  precipitating  parti¬ 
cles. 

The  SSJ/4  detectors  on  each  satellite  consist  of  an  array  of 
four  electrostatic  analyzers  that  together  measure  electrons 
and  ions  over  an  energy  range  of  30  eV  to  30  keV  in  20 
channels  logarithmically  spaced  in  energy.  A  20-point  spec¬ 
trum  is  returned  each  second  for  both  electrons  and  ions.  A 
complete  description  of  the  SSJ/4  instrument  is  given  by 
Hardy  et  al.  [1984].  The  instruments  were  extensively  cali¬ 
brated  prior  to  flight  and  cross-calibrated  in  flight. 

Although  both  DMSP  satellite  orbits  are  sun  synchronous, 
each  day  they  sweep  out  a  wide  range  in  corrected  geomag¬ 
netic  latitude  (MLAT)  and  magnetic  local  time  (MLT)  coor¬ 
dinates  owing  to  the  offset  between  geographic  and  geomag¬ 
netic  poles  in  combination  with  the  Earth’s  rotation.  Figure 
1  shows  the  data  coverage  (combined  hemispheres)  in 
MLAT-MLT  coordinates  resulting  from  hundreds  of  DMSP 
F6  and  F7  orbits.  MLAT  is  marked  in  1(T  increments  from 
50°  to  90°  by  the  concentric  circles.  The  dots  along  the  50° 
MLAT  circle  appear  evf  0.5  hour  in  MLT.  Coverage  of  the 
auroral  region  is  nearly  complete,  with  a  data  gap  at  low 
magnetic  latitudes  in  the  post.ioon  and  postmidnight  local 
time  sectors. 
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Fig.  I.  The  data  coverage  from  DMSP  F6  and  F7  satellites  for 
combined  hemispheres  displayed  in  a  corrected  geomagnetic  lati¬ 
tude-magnetic  local  time  coordinate  system. 


3.  Data  Binning  Procedure 

We  have  constructed  identically  structured  but  separate 
data  bases  for  both  electrons  and  ions.  Our  procedure  f  • 
binning  and  averaging  the  data  is  identical  to  that  of  papers 
I  and  2  except  that  instead  of  sorting  the  data  by  the  KP 
index,  we  sort  it  by  paired  values  of  IMF  B .  (in  geocentric 
solar  magnetospheric  coordinates)  and  VbW.  The  matrix  of 
sorting  parameters  has  six  divisions  in  B,  ranging  from  - 10 
to  + 10  nT  and  five  divisions  in  Vsw  ranging  from  200  to  800 
km/s,  resulting  in  30  paired  ranges  of  interplanetary  param¬ 
eters.  The  bin  intervals  for#,  and  are  listed  in  Table  1. 
Particle  spectra  sorted  into  each  of  the  30  categories  are 
further  ordered  into  a  spatial  grid  (map)  in  corrected  geo¬ 
magnetic  latitude  (MLAT)  and  local  time  (MLT).  The  reso¬ 
lution  of  the  map  is  2°  for  50°-60°  and  80°-90°  MLAT  and  1° 
for  60°-80°  MLAT;  it  is  0.5  hour  in  MLT.  The  magnetic 
coordinates  are  determined  by  mapping  the  satellite  coordi¬ 
nates  down  to  100  km  using  the  International  Geophysical 
Reference  Field  (1982)  and  the  Jensen-Cain  magnetic  field 
models. 

Previous  work  [Baker  et  al.,  1983,  and  references  therein] 
has  shown  that  interplanetary  parameters  are  most  highly 
correlated  with  magnetospheric  activity  parameters  when  a 
time  delay  is  assumed.  This  time  delay  ranges  from  0.25  to  1 


TABLE  1.  Bin  Intervals 


Bz,  nT 

v 

V  SW* 

km/s 

Interval 

Average 

Interval 

Average 

- 10.0 —  3.0 

-4.8 

200-375 

346 

-3.0— 1.5 

-2.2 

375-450 

408 

-1. 5-0.0 

-0.7 

450-525 

485 

0.0-4- 1.5 

+0.7 

525-625 

572 

+ 1 .5“+ 3.0 

+  2.1 

623-800 

677 

+3.0-+ 10.0 

+4.5 

hour  depending  on  the  level  of  activity  in  the  interplanetary 
medium  and  can  be  explained  in  terms  of  the  propagation 
time  of  solar  wind  perturbations  as  well  as  the  finite  response 
time  of  the  magnetosphere.  A  study  we  performed  using 
1826  hours  of  data  showed  that  ordering  1-hour  averages  of 
the  AE  index  by  B.  and  VbW  gives  less  scatter  for  1-hour  than 
for  0-hour  lag  time.  We  therefore  chose  a  1-hour  time  lag 
between  the  hourly  averaged  interplanetary  data  obtained 
from  the  National  Space  Science  Data  Center  (NSSDC)  data 
base  [Couzens  and  King,  1986,  1989]  and  the  precipitating 
particle  data.  Because  the  particle  data  resolution  is  1  s,  the 
stated  time  lag  of  1  hour  is  used  in  an  average  sense  only.  For 
each  hour  a  pair  of  IMF  B,  and  VbW  values  exist;  all  1-s 
particle  spectra  in  the  following  1-hour  interval  are  sorted 
according  to  that  pair  of  hourly  averaged  interplanetary 
parameters. 

Simultaneous  IMF  and  plasma  data  were  available  for 
approximately  two  thirds  of  time  between  January  1983  and 
July  1987,  resulting  in  a  total  of  13,263  pairs  of  hourly  values. 
The  data  bases  for  electrons  and  ions  each  contain  a  total  of 
34  million  spectra  distributed  between  the  30  individual 
maps,  each  map  consisting  of  between  0.2  and  2.7  million 
spectra.  The  average  spectrum  in  each  spatial  bin  is  derived 
from  hundreds  to  thousands  of  1-s  spectra  sampled  from  tens 
to  hundreds  of  individual  passes.  Within  a  given  bin  the 
independence  between  individual  1-s  spectra  varies  widely 
with  their  temporal  proximity  to  one  another.  The  large 
number  of  separate  passes  contributing  to  each  bin  mini¬ 
mizes  the  impact  that  the  autocorrelation  between  tempo¬ 
rally  clustered  1-s  spectra  rnay  have  on  the  average  spec¬ 
trum. 

Table  2  summarizes  the  bin  statistics  and  parameteriza¬ 
tion.  Included  for  each  of  the  30  separations  are  the  total 
number  of  spectra  used  in  each  map  and  the  number  of 
paired  hourly  values  of  B.  and  V'bW  along  with  their  averages 
and  standard  deviations  (a).  The  table  is  arranged  so  that  the 
six  B.  bins  are  listed  repeatedly  for  increasing  values  of  fsw. 
To  simplify  the  bookkeeping,  a  single  weighted  average 
value  of  each  of  the  B.  bins  was  used  independently  of  Vsw. 
Likewise,  a  single  value  for  each  of  the  Vsw  bins  was  used 
independently  of  B..  These  average  bin  values  are  listed  in 
Table  1 ,  along  with  their  corresponding  bin  intervals.  The 
differences  between  these  averages  and  the  averages  listed  in 
Table  2  are  minor  and  inconsequential  to  the  final  results. 
Also  included  in  Table  2  are  the  average  values  of  Kp  and 
the  solar  wind  dynamic  pressure  associated  with  each  map. 
The  average  Kp  values  will  be  used  in  a  later  section  for 
comparing  the  current  work  to  our  prevbus  Kp  statistical 
studies  (papers  1  and  2).  The  average  dynamic  pressure 
(proton  mass  times  solar  wind  density  times  Vsvv  squared) 
has  been  determined  from  the  NSSDC  data  base  and  will  be 
referred  to  in  the  discussion  (section  6). 

Two  points  should  be  made  concerning  the  map  separa¬ 
tions.  One  is  that  the  extreme  B.  negative  and  positive  bins 
(7  nT  wide)  are  much  larger  than  the  four  moderate  B.  bins 
(1.5  nT  wide).  These  wide  bins  were  necessitated  by  the 
paucity  of  strong  B.  conditions  and  the  requirement  for 
adequate  bin  statistics.  For  these  wide  bins  the  resulting 
standard  deviations  in  B.  are  about  3  times  those  found  in 
the  four  moderate  B.  bins  (see  Table  2).  The  second  point  is 
that  there  is  an  apparent  minimum  number  of  spectra  re¬ 
quired  to  produce  a  smoothly  varying  and  statistically  sig¬ 
nificant  map.  We  find  this  number  to  be  between  0.25  and  0.5 
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TABLE  2.  Bin  Averages  and  Statistics 


Map 

Number 

Spectra, 
x  106 

Hours 

Bz,  nT 

V5W,  km/s 

Kp 

Pk.  keV/cm3 

Average 

Sigma 

Average 

Sigma 

Average 

Sigma 

Average 

Sigma 

0.81 

290 

-4.5 

1.3 

346. 

20. 

2.8 

0.8 

21.6 

10.1 

2 

1.20 

405 

-2.2 

0.4 

346. 

20. 

1.9 

0.8 

19.6 

9.0 

3 

2.25 

780 

-0.7 

0.4 

345. 

20. 

1.4 

0.8 

18.5 

8.3 

4 

2.25 

811 

0.7 

0.4 

345. 

21. 

1.1 

0.8 

18.1 

8.2 

5 

1.40 

501 

2.1 

0.4 

347. 

21. 

1.1 

0.8 

17.8 

8.1. 

6 

1.02 

360 

4.5 

1.4 

349. 

22. 

1.3 

0.9 

20.4 

9.1 

7 

1.12 

441 

-4.9 

1.5 

408. 

21. 

3.6 

1.1 

22.1 

11.9 

8 

1.45 

532 

-2.2 

0.4 

409. 

21. 

2.7 

0.9 

19.0 

10.1 

9 

1.60 

923 

-0.8 

0.4 

408. 

21. 

2.0 

0.9 

17.9 

10.7 

10 

2.66 

1,007 

0.7 

0.4 

409. 

21. 

1.7 

0.9 

17.1 

9.6 

11 

1.92 

697 

2.1 

0.4 

408. 

20. 

1.6 

0.9 

17.9 

9.1 

12 

1.52 

562 

4.6 

1.5 

407. 

21. 

1.8 

1.1 

22.5 

13.9 

13 

0.72 

311 

-5.0 

1.7 

485. 

22. 

4.3 

1.1 

26.9 

16.1 

14 

0.80 

332 

-2.2 

0.4 

485. 

22. 

3.2 

0.8 

19.3 

10.9 

15 

1.38 

553 

-0.7 

0.4 

484. 

23. 

2.7 

0.9 

17.7 

8.5 

16 

1.65 

653 

0.7 

0.4 

484. 

22. 

2.2 

1.1 

17.3 

10.8 

17 

1.13 

421 

2.1 

0.4 

487. 

22. 

2.0 

1.2 

17.5 

11.2 

18 

0.86 

348 

4.7 

1.7 

486. 

20. 

2.4 

1.2 

25.8 

14.3 

19* 

0.42 

197 

-4.7 

1.5 

569. 

25. 

4.6 

1.0 

27.7 

13.3 

20 

0.61 

255 

-2.3 

0.5 

573. 

29. 

3.9 

0.9 

21.3 

8.6 

21 

1.03 

427 

-0.7 

0.4 

572. 

28. 

3.4 

1.1 

20.5 

9.2 

22 

1.23 

502 

0.7 

0.4 

573. 

28. 

2.9 

1.0 

19.0 

8.2 

23 

0.98 

410 

2.1 

0.4 

571. 

28. 

2.6 

1.1 

1.9.3 

8.3 

24 

0.65 

284 

4.4 

1.4 

572. 

29. 

3.0 

1.2 

26.2 

13.4 

25* 

0.16 

71 

-4.4 

1.4 

668. 

36. 

5.0 

1.1 

25.9 

11.1 

26» 

0.35 

149 

-2.2 

0.4 

675. 

40. 

4.4 

0.9 

23.3 

10.3 

27 

0.66 

256 

-0.8 

0.4 

679. 

44. 

3.9 

0.9 

22.1 

7.7 

28 

0.87 

368 

0.7 

0.4 

679. 

42. 

3.5 

0.9 

21.1 

8.0 

29 

0.63 

261 

2.2 

0.4 

675. 

38. 

3.2 

1.0 

21.8 

10.3 

30* 

0.42 

156 

4.3 

1.3 

677. 

36. 

3.0 

1.2 

27.3 

15.9 

33.75t 

13,2631" 

•Maps  omitted  from  study. 
tTotals. 


million  spectra,  the  exact  level  depending  on  the  intrinsic 
variability  of  the  data,  which  in  turn  depends  on  the  activity 
and  the  width  of  the  bins.  We  note  that  maps  19,  25,  26,  and 
30  consist  of  fewer  than  0.5  million  spectra  each,  with  the 
majority  of  their  average  spectra  defined  by  less  than  10 
passes  of  data  and  representing  an  insignificant  number  of 
truly  independent  1-s  spectra.  Of  these  maps,  all  but  number 
26  correspond  to  a  wide  Bz  bin  which  compounds  the 
difficulty  of  providing  meaningful  average  hemispheric 
fluxes.  Therefore  aside  from  including  the  data  from  these 
four  maps  in  Tables  2  and  3,  they  have  been  omitted  from  the 
study. 

4.  Definition  of  Average  Spectral  Moments 

Average  spectra  are  determined  from  the  1-s  spectra 
accumulated  in  each  spatial  bin  for  each  map.  From  these 
average  spectra  (combined  north  and  south  poles),  we  de¬ 
termine  for  both  ions  and  electrons  the  total  precipitating 
number  and  energy  flux  into  the  hemisphere  each  second 
urder  different  average  interplanetary  conditions.  To  pro¬ 
vide  the  best  estimate  of  the  total  hemispheric  particle  flux 
input  in  the  statistical  results  presented  here,  we  fill  the  data 
gaps  indicated  in  Figure  I  by  interpolation.  Upon  examina¬ 


tion  of  the  flux  maps  created  in  this  study  (not  shown),  it  is 
found  that  flux  levels  adjacent  to  the  data  gap  in  the 
postmidnight  sector  are  above  backround  for  s58°  MLAT. 
Al  these  latitudes,  interpolation  between  real  flux  data  is 
limited  to  s2.5  hours  in  MLT.  Fluxes  adjacent  to  the 
postnoon  sector  data  gap  are  at  backround  level,  and  inter¬ 
polation  across  the  gap  is  of  no  consequence  to  the  final 
results. 

In  about  30%  of  the  average  ion  (electron)  spectra,  partic¬ 
ularly  those  determined  for  high  levels  of  geomagnetic 
activity  on  the  evening  (morning)  side  of  the  oval,  a  substan¬ 
tial  portion  of  the  high-energy  tail  of  the  spectra  lies  beyond 
the  30-kcV  upper  bound  of  the  SSJ/4  instrument.  This  can 
underestimate  (more  significantly  for  ions  than  electrons)  the 
total  precipitating  integral  number  and  energy  flux.  To 
correct  for  this,  we  selectively  extrapolate  average  spectra 
to  three  additional  logarithmically  spaced  energies  (45, 65, 90 
keV)  assuming  the  tails  of  the  spectra  to  be  Maxwellian 
extensions  of  the  high-energy  measuremei's  (see  paper  2  for 
a  detailed  discussion  of  this  conseivative  correction). 

From  the  average  spectra  we  calculate  the  integra*  energy 
flux  (JEflux  'n  keV/s  cm2  sr)  and  integral  number  flux 
(Afflux  in  particles/s  cm2  sr)  using  the  equations 


Brautigam  et  al.;  Effect  of  IMF  Bz  and  Vsw  on  Precipitation 


5529 


^FLUX  -  E\  ‘j(E  l)  *  (^2  -El) 


22 

+  2  ErJlEd- 


4-2 


(Ei  +  ;  Ei  -  i ) 
2 


+  £23  'JiEii)  •  (£’23  “  ^22)  (1) 

(is  —  £■  ) 

Afflux  =y(£i)  •  (£2  -  £|)  +  2  J(£.)  • 

i  =  2  1 

+  j(E2i)-(E2*-E22)  (2) 

where  „>(£,)  is  the  differential  ion  or  electro'  .nber  flux  in 
units  of  particles/(cm2  s  sr  keV)  for  the  n  energy  channel 
and  Ej  is  the  central  energy  ir.  units  of  keV.  The  summation 
runs  over  the  20  measured  channels  plus  the  three  extrapo¬ 
lated  channels.  We  integrate  NFlux  afid  £flux  over  each 
spatial  map  (i.e.,  above  50°  MLAT)  to  obtain  for  each  of  the 
30  average  interplanetary  conditions  the  precipitating  parti¬ 
cle  hemispheric  integral  number  flux  (# nflux.  in  particles/(s 
sr))  and  hemispheric  integral  energy  flux  (#£tLUX>  in  giga¬ 
watts).  As  in  paper  2,  we  assume  that  the  precipitating 
particle  distributions  are  isotropic  and  thus  introduce  a 
factor  of  pi  in  converting  HEfLux  to  gigawatts.  Throughout 
the  text  we  use  the  convention  of  preceding  Nflux  and 
£flux  wit*1  an  "H"  when  referring  to  the  spatially  inte¬ 
grated  hemispheric  flux. 

Approximately  30%  of  the  average  ion  and  electron  spec¬ 
tra  were  extrapolated.  In  the  majority  of  the  extrapolated  ion 
spectra,  the  £r.'LUX  was  increased  by  less  than  50%,  while 
the  Afflux  increased  by  less  than  10%.  The  effect  of 
extrapolation  on  the  electron  spectra  is  less  significant  with 
the  £Flux  increasing  by  less  than  20%  and  the  WFlux 
increasing  by  less  than  10%.  These  increases  in  integral 
fluxes  for  individual  spectra  translate  to  a  5-90%  (0-9%) 
increase  in  the  ion  (electron)  hemispheric  energy  flux  and  a 
1-12%  (0-0.4%)  increase  in  the  ion  (electron)  hemispheric 
number  flux,  depending  on  activity  level.  In  listing  hemi- 
spherically  integrated  values,  both  those  from  unextrapo¬ 
lated  and  extrapolated  spectra  are  given. 

For  additional  details  on  data  reduction  techniques  (cross¬ 
satellite  normalization,  smoothing,  and  extrapolation),  we 
refer  the  reader  to  papers  1  and  2. 

5.  Results 

The  spatial  patterns  in  the  £Flux  and  NFlux  maF! 
generated  in  this  study  are  similar  to  those  shown  previously 
in  papers  1  and  2;  only  the  quantitative  details  differ.  We 
therefore  refer  the  reader  to  the  previous  papers  for  a 
discussion  of  the  morphological  features  of  auroral  precipi¬ 
tation,  and  in  this  paper  we  concentrate  on  the  total  hemi¬ 
spheric  energy  and  number  flux  for  the  precipitating  ions  and 
electrons  and  their  dependence  on  the  solar  wind  speed  and 
the  IMF  Bz  component.  Independent  of  the  spatial  distribu¬ 
tion  over  the  auroral  zone  of  the  energy  and  number  flux,  we 
take  the  integrated  hemispheric  inputs  ( HEnm  and  HNnux) 
to  be  a  measure  of  the  character  and  strength  of  the  solar 
wind-magnetosphere  coupling.  Detailed  analysis  of  the  de¬ 
pendence  of  the  precipitation  patterns  and  the  average 
spectra  on  Vsw  and  IMF  Bz,  especially  as  they  pertain  to  the 


various  auroral  regions  such  as  the  boundary  layers  and  the 
central  plasma  sheet,  will  be  addressed  in  a  future  paper. 

5.1.  Hemispheric  Flux  as  a  Function  of  Bz  and  Vsw 

The  and  as  determined  by  integrating  the 

extrapolated  average  spectra  over  the  spherical  area  of  the 
auroral  region,  are  listed  in  Table  3  and  plotted  in  Figures 
2-5.  The  fluxes  determined  from  the  unextrapolated  spectra 
are  also  included  in  parentheses  in  Table  3.  In  Figures  2-5 
the  electron  results  are  plotted  on  the  left,  and  the  ion  results 
are  shown  on  the  right.  All  scales  are  linear.  The  hemi¬ 
spheric  number  flux  plots  have  different  multiplicative  fac¬ 
tors  for  electrons  (1025)  and  ions  ( 1 0 23 ) .  We  have  plotted  the 
data  in  two  ways  to  illustrate  the  two  principal  characteris¬ 
tics  of  the  B,  and  Vsw  variation.  //W(u.x  and  //£|LUX  are 
plotted  against  B.  for  fixed  Vsw  in  Figures  2  and  3,  respec¬ 
tively,  and  are  plotted  against  Vsw  for  fixed  Bz  in  Figures  4 
and  5,  respectively.  Different  symbols,  as  defined  in  the 
symbol  key,  are  used  for  each  different  Bz  and  Vsw  case.  As 
was  previously  mentioned  in  section  3,  values  from  four  of 
the  initial  30  maps  have  been  omitted  and  will  not  be  used  in 
any  of  the  following  analyses  owing  to  their  poor  statistics. 

The  electron  and  ion  hemispheric  number  flux  (Figure  2) 
and  ion  hemispheric  energy  flux  (Figure  3)  vary  similarly 
with  B , :  decreasing  from  an  absolute  maximum  for  strong  Bz 
negative  to  a  minimum  for  weak  to  moderate  Bz  positive, 
followed  by  an  increase  to  a  second  maximum  for  strong  Bz 
positive.  The  electron  hemispheric  energy  flux  (Figure  3) 
falls  off  rapidly  from  its  maximum  for  strong  B.  negative  to 
weak  Bz  positive,  gradually  leveling  off  to  its  minimum  for 
strong  B .  positive.  The  electron  and  ion  energy  fluxes  show 
the  largest  overall  increase  from  their  minimum  to  their 
maximum  for  strong  Bz  negative,  by  a  factor  of  5.8  and  3.7, 
respectively.  This  occurs  for  the  lowest  V5W  bin,  with  the 
factors  decreasing  for  higher  Vsw.  The  electron  and  ion 
number  fluxes  show  smaller  rises,  both  increasing  by  a  factor 
of  approximately  2.5  and  remaining  near  this  value  for  all 
Vsw.  For  Bz  positive  of  the  minimum,  all  but  the  electron 
HElLm  increase.  The  increase  is  largest  for  the  ions  with  an 
average  rise  of  a  factor  of  1.41  and  1.42  for  the  ion  HEflvx 
and  HNtwx,  respectively.  The  electron  H W|XUX  increases  on 
average  by  a  factor  of  1.19. 

Figures  4  and  5  illustrate  that  for  each  B.  value  the  four 
hem^pherically  integrated  fluxes  exhibit  a  near-lmear  vari¬ 
ation  with  Vsw.  Again  the  electron  and  ion  HEyLVX  show  the 
largest  variation,  increasing  by  an  average  factor  of  4.08  and 
2.11,  respectively,  over  the  full  velocity  range.  The  electron 
and  ion  HNtlm  increase  by  a  factor  of  1.82  and  1.16, 
respectively. 

The  similarity  in  the  variation  of  the  four  quantities  with 
B,  and  Vsw  suggests  that  the  variation  can  be  fit  by  either  of 
two  relatively  simple  functional  forms.  In  each  case,  the  B. 
functionality  is  modeled  by  a  parabolic  least  squares  fit  to  the 
data,  and  the  Vsw  functionality,  by  a  linear  least  squares  fit. 
In  the  first  case,  the  two  functions  are  summed: 

f(Bz,  Vsw)  =  a(Bz  -  h )2  +  [cVsw  +  d]  (3) 
while  in  the  second  case,  they  are  multiplied: 

f(B;,  =  [a(Bz  -  h )2  +  l][rVw  +  d] |.  (4) 

Although  equations  (3)  and  (4)  have  four  free  parameters 
each,  the  functional  forms  differ.  Equation  (3)  is  the  simpler, 
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TABLE  3.  Hemispheric  Fluxes  (Electron  and  Ion) 


Map 

Number 

BrnT 

Vsw,  km/s 

Electron 

tfW’GW 

Electron 

x,o"%UsXr’r' 

Ion 

GW 

Ion 

Kdum' 
x  1023  (s  srl  ■ 

1 

-4.5 

346. 

35.00  (34.50) 

3.83  (3.83) 

3.71  (2.63) 

7.03  (6.61) 

2 

-2.2 

346. 

20.80  (20.50) 

2.68  (2.68) 

2.19(1.77) 

4.78  (4.61) 

3 

-0.7 

346. 

11.80  (11.70) 

2.01  (2.01) 

1.37  (1.20) 

3.65  (3.58) 

4 

0.7 

346. 

7.18(7.17) 

1.74  (1.74) 

1.18(1.05) 

3.65  (3.60) 

5 

2.2 

346. 

5.71  (5.70) 

1.65  (1.65) 

1.02(0.96) 

3.60  (3.58) 

6 

4.5 

346. 

6.23(6.17) 

2.01  (2.01) 

1.38  (1.28) 

5.44  (5.39) 

7 

-4.5 

408. 

49.50  (48.50) 

5.05  (5.04) 

5.05  (3.16) 

8.10  (7.39) 

8 

-2.2 

408. 

26.80  (26.40) 

3.03  (3.03) 

2.56  (1.83) 

4.88  (4.601 

9 

-0.7 

408. 

18.60(18.40) 

2.55  (2.55) 

1.93  (1.49) 

4.01  (3.84) 

10 

0.7 

408. 

11.50(11.40) 

2.04  (2.04) 

1.56  (1.23) 

3.55  (3.42) 

11 

2.2 

408. 

9.58  (9.54) 

2.08  (2.08) 

1.36(1.19) 

3.82  (3.75) 

12 

4.5 

408. 

8.82  (8.79) 

2.41  (2a;) 

1.78  (1.55) 

5.35  (5.26) 

13 

-4.5 

485. 

65.90  (63.30) 

5.87  (5.86) 

5.72  (3.21) 

8.62  (7.71) 

14 

-2.2 

485. 

37.60  (36.50) 

3.54  (3.54) 

2.68  (1.76) 

4.83  (4.48) 

Ic 

-0.7 

485. 

27.40  (26.90) 

2.87  (2.87) 

2.09(1.45) 

3.87  (3.62) 

16 

0.7 

485. 

15.90  (15.80) 

2.30  (2.30) 

1.86(1.30) 

3.61  (3.40) 

17 

2.2 

485. 

12.80(12.70) 

2.23  (2.23) 

1.52  (1.23) 

3.60  (3.48) 

18 

4.5 

485. 

11.60(11.50) 

2.77  (2.77) 

2.21  (1.79) 

5.46  (5.30) 

19* 

-4.5 

572. 

68.40  (64.60) 

5.93  (5.92) 

5.70  (3.05) 

8.53  (7.58) 

20 

-2.2 

572. 

46.90  (44.70) 

4.07  (4.06) 

3.31  (1.98) 

5.74  (5.25) 

21 

-0.7 

572. 

38.80  (37.40) 

3.54  (3.54) 

2.83  (1.65) 

4.64(4.21) 

22 

0.7 

572. 

25.40  (24.70) 

2.87  (2.86) 

2.30(1.41) 

3.78  (3.45) 

23 

2.2 

572. 

18.70  (18.30) 

2.77  (2.77) 

1.89  (1.32) 

3.60  (3.38) 

24 

4.5 

572. 

18.00(17.80) 

3.18(3.18) 

2.85  (1.95) 

5.48(5.15) 

25* 

-4.5 

677. 

65.50  (60.00) 

5.42  (5.40) 

4.61  (2.27) 

6.64  (5.86) 

26* 

-2.2 

677. 

60.20  (55.20) 

4.66  (4.64) 

3.28  (1.89) 

5.64  (5.13) 

27 

-0.7 

677. 

43.50(41.10) 

3.82  (3.81) 

2.86(1.64) 

4.66  (4.22) 

28 

0.7 

677 

36.30  (34.70) 

3.49  (3.48) 

2.42(1.42) 

3.99  (3.63) 

29 

2.2 

677. 

25.00  (24.30) 

3.11  (3.11) 

2.36  (1.43) 

4.04  (3.70) 

30* 

4.5 

677. 

20.30  (19.80) 

3.13(3.13) 

2.76  (1.77) 

4.85  (4.49) 

Unextrapolated  values  are  listed  in  parentheses. 
*Maps  omitted  from  study. 


containing  only  the  terms  B\,  Bz,  and  ViW.  Equation  (4) 
contains  the  additional  cross  terms,  VSVlBz  and  VhWB:.  The 
four  constants,  a  through  d,  for  the  two  methods  of  fitting 
are  listed  in  Table  4.  Because  there  are  four  free  parameters 


in  each  method,  the  quality  of  fit  should  be  approximately 
the  same.  A  quantitative  comparison  between  the  hemi¬ 
spheric  fluxes  listed  in  Table  3  and  those  calculated  by 
applying  the  fit  coefficients  in  Table  4  (methods  1  and  2)  to 
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Fig.  2.  The  hemispheric  number  flux  (///vin,x).  in  particles/(s  sr).  plotted  as  a  function  of  the  bin-averaged  value  of 

IMF  lL.  (Left)  Electrons.  (Right)  Ions. 
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Fig.  3.  The  hemispheric  energy  flux  (//£|  J  UJ,  in  gigawatts,  plotted  as  a  function  of  the  bin-averaged  value  of  IMF  B. 

(Left)  Electrons.  (Right)  Ions. 


equations  (3)  and  (4),  respectively,  show  the  following: 
Method  1  results  in  71%  (94%)  of  the  hemispheric  fluxes 
being  fit  to  within  10%  (20%)  of  the  data.  The  correlation 
coefficient  for  a  linear  regression  run  between  the  data  and 
fitted  points  for  each  category  of  fluxes  (Hehux  and  HNum 
for  both  electrons  and  ions)  averaged  0.95.  Method  2  yielded 
slightly  better  fits,  with  80%  (99%)  of  fitted  points  within  10% 
(20%)  of  the  data.  The  average  correlation  coefficient  for 
method  2  was  0.98.  A  more  detailed  discussion  of  the  fitting 
methods  and  results  are  provided  as  an  appendix  to  this 
report. 

5.2.  Comparison  of  Ion  and  Electron 
Hemispherically  Integrated  Fluxes 

A  cursory  examination  of  Table  3  shows  that  the  ion 
he„vx  Wnilvx)  >s  always  less  than  the  electron  //£hux 
Mv„ux)  by  about  1  (2)  orders  of  magnitude.  However,  a 


closer  examination  of  the  ratios  between  ion  and  electron 
hemispheric  fluxes  (#£nux  and  Hnuux)  show  important 
variations  with  respect  to  B,  and  Vsw.  Figure  6  left  (right) 
shows  the  ratio  of  ion  to  electron  (//^n  ux)  plotted  as 

a  function  of  Bz  for  fixed  Vsw. 

The  largest  ratio  (0.22)  of  ion  to  electron  HEnm  (e.g., 
when  the  two  are  most  comparable)  occurs  for  the  lowest 
velocities  and  for  Bz  strongly  northward.  The  ratio  de¬ 
creases  sharply  as  B ,  decreases  toward  zero  and  levels  off  to 
an  approximately  constant  value  for  Bz  negative.  This  trend 
occurs  at  each  velocity  level.  The  ratios  vary  in  total  by  a 
factor  of  4  from  0.054  to  0.223.  For  a  fixed  B.  the  ratio 
decreases  with  increasing  velocity  in  a  roughly  linear  man¬ 
ner.  The  rate  of  decrease  is  greater  for  B .  positive  than  for 
B,  negative. 

A  comparable  variation  is  seen  in  the  Hnhvx  ratio-  The 
largest  values  are  seen  for  the  lowest  velocity  and  the  largest 
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Fig.  4.  The  hemispheric  number  flux  (/fjvILtiX).  in  particles/(s  sr),  plotted  as  a  function  of  the  bin-averaged  value  of 

Vsw.  (Left)  Electrons.  (Right)  Ions. 
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Fig.  5.  The  hemispheric  energy  flux  in  gigawatts,  plotted  as  a  function  of  the  bin-averaged  value  of  Vsw 

(Left)  Electrons.  (Right)  Ions. 


Bz  positive.  For  B.  positive  the  value  decreases  rapidly  for 
decreasing  B,  for  all  Vsv/.  For  Bz  negative  the  value  is  approx¬ 
imately  constant  for  a  given  Vsw.  Again  there  is  a  larger 
variation  with  velocity  for  Bz  positive  than  for  B.  negative. 

5.3.  Variation  of  Kp  With  B.  and  Vsw 
As  a  point  of  comparison  to  our  previous  studies,  we  have 
used  3-hour  Kp  averages  to  determine  the  variation  in  Kp 
with  B.  and  Vs>v  for  the  data  set  used  here.  Similar  studies  of 
magnetic  activity  variation  with  solar  wind  parameters  have 
been  performed,  and  those  that  most  closely  parallel  this 
analysis  are  given  by  Berthelier  [1976]  and  Maezawa  and 
Murayama  [1986].  The  results  are  tabulated  in  Table  2  and 
plotted  in  Figure  7,  where  the  separations  in  B.  and  Vsw  are 
the  same  as  for  the  HEnux  and  HNrlm  studies.  One  notes 
that  Kp  exhibits  the  same  general  variation  with  B.  and  Vsw 
as  do  HEiu)X  and  Hnuvx-  Kp  has  a  minimum  value  for  B. 
weak  to  moderately  positive  and  increases  slightly  for  Bz 
strongly  positive.  Kp  increases  approximately  linearly  with 
V5* .  This  similarity  in  behavior  is  not  surprising,  since  Kp 
should  respond  to  the  general  level  of  geomagnetic  activity 


in  the  oval,  which  in  turn  should  be  proportional  to  the  total 
number  and  energy  flux  input  into  the  oval. 

Figure  7  also  shows  that  any  separation  of  the  data  by  Kp 
will  mix  periods  of  different  polarity  in  Bz  and  have  a  broad 
range  in  Vsw.  In  Figure  8  we  show  the  distribution  of  B,  for 
Kp  =  0,  2,  4,  and  6,  over  the  same  time  period  from  which 
our  data  base  is  made.  Plotted,  for  four  Kp  levels,  is  the 
fraction  of  time  that  Bz  falls  within  each  of  1.5-nT-wide  bins 
over  the  range  - 10.5  to  + 10.5  nT.  The  figure  illustrates  that 
even  for  low  activity  (Kp  =  0),  the  IMF  B,  component  is 
negative  12%  of  the  time.  Conversely,  for  very  active 
periods  (Kp  =  6),  the  IMF  Bz  component  is  positive  23%  of 
the  time.  For  more  moderate  periods  (Kp  =  2  and  Kp  =  4) 
the  percentage  of  opposite  B.  polarities  is  approximately 
equal.  Thus  a  separation  by  Kp  poorly  reflects  groupings  of 
those  solar  wind  parameters  that  order  the  degree  of  the 
solar  wind-magnetospheric  coupling. 

6.  Discussion 

We  have  shown  that  the  average  hemispheric  input  of 
particle  energy  and  particle  flux  for  precipitating  electrons 


TABLE  4. 

Hemispheric  Flux  Fit  Coefficients 

Electrons 

Ions 

Coefficient* 

Energy  Flux 

Number  Flux 

Energy  Flux 

Number  Flux 

a,  x/nT2 

0.59 

Method  1 

0.06 

0.07 

0.14 

b,  nT 

4.45 

2.33 

2.00 

0.96 

c,  x/(km/s) 

0.092 

0.005 

0.004 

0.002 

d,  x 

-34.20 

-0.31 

-0.38 

2.66 

-die,  km/s 

372. 

62. 

95. 

-1330. 

a,  1/nT 

0.08 

Method  2 

0.03 

0.06 

0.04 

b,  nT 

3.31 

2.00 

1.88 

0.95 

c,  .v/(km/s) 

0.041 

0.004 

0.004 

0.001 

X 

-8.65 

-0.21 

-0.32 

3.29 

-die,  km/s 

211. 

-53. 

80. 

-3290. 

*Variable  ,x  is  in  gigawatts  (energy  flux)  or  (s  sr)"1  (number  flux). 
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Fig.  6.  The  ratio  of  ion  to  electron  hemispheric  flux  (left,  //£flux;  right,  Usri  ux)  plotted  as  a  function  of  Bz. 


ar, 1  ions  in  the  auroral  region  varies  in  a  consistent  and 
relatively  simple  manner  with  the  magnitude  of  the  IMF  B. 
component  and  with  the  solar  wind  velocity  Vsw.  In  this 
section  we  examine  the  implications  of  these  variations  in 
more  detail.  Specifically,  we  address  the  following  issues:  (1) 
the  statistical  significance  of  the  increase  in  //E[L(jx  and 
HNfLVX  for  large  B,  positive,  (2)  the  possible  role  of  solar 
wind  pressure  in  the  increase  for  Bz  posi'ive,  (3)  the  signif¬ 
icance  of  the  functional  form  that  results  from  fitting  the 
data,  (4)  the  role  of  the  solar  wind  velocity  in  controlling  the 
energy  transfer,  and  (5)  the  comparison  of  these  results  to 
earlier  studies. 

6. 1 .  Possible  Bias  in  the  Bz  Distribution 

We  first  address  whether  the  increase  in  HEftvx  and 
HNflvx  (predominantly  for  ions)  for  increasing  Bz  positive  is 
a  residual  effect  of  preceding  Bz  negative  intervals  or  is  a  real 
effect  associated  with  strong  Bz  positive.  Nakai  and  Kamide 


[1983]  have  shown  that  the  magnetospheric  response  to  Bz 
positive  as  measured  by  the  motion  of  the  equatorward 
boundary  of  the  oval  is  much  slower  than  for  Bz  negative. 
This  implies  that  for  several  hours  after  a  northward  turning 
of  Bz ,  the  //£mjx  and  HNfLm  could  be  affected  in  part  by  the 
residual  influence  of  the  increasingly  more  remote  Bz  nega¬ 
tive  time  period.  The  only  way  this  could  explain  the 
observed  rise  in  HEfwx  and  HNflw  with  increasing  B. 
positive  would  be  if  the  probability  that  a  period  of  large 
negative  Bz  prior  to  a  period  of  positive  Bz  increases  for 
increasing  magnitude  of  positive  Bz. 

To  investigate  this,  we  used  the  IMF  data  covering  the 
time  period  of  1983-1986.  From  this  set  we  extracted  those 
hourly  Bz  positive  values  which  belonged  to  a  time  interval 
that  had  continuous  coverage  running  backward  in  time  to  an 
hour  where  Bz  fell  below  -2  nT.  These  cases  were  binned 
according  to  their  magnitude  (two  bins  defined  by  0  to  +3 
nT,  +3  to  + 10  nT),  as  well  as  by  the  number  of  hours  (up  to 


-5-3-1  I  3  5  300  400  500  600  700 
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Fig.  7.  The  average  value  of  Kp  for  each  of  the  26  bins  (Table  2)  plotted  as  a  function  of  (left)  bin-averaged  B  ■  values 

and  (right)  bin-averaged  Vsw  values. 
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Fig.  8.  Distribution  of  the  fraction  of  time  that  B  falls  within 
each  of  the  various  B.  bins  (or  Kp  =  0,  0+ ;  2-,  2,  2+.:  4-.  4,  4+; 
and  a6-. 


13)  between  them  and  the  nearest  preceding  negative  B.. 
This  sort  resulted  in  4032  (1418)  hours  in  the  weak  (strong) 
B ,  positive  bin.  The  results  of  this  analysis  are  plotted  in  Figure 
9.  Plotted  as  a  function  of  the  hours  preceding  a  given  B. 
positive  measurement  for  both  B.  bins  is  the  fraction  of  cases 
where  a  Bz  of  less  than  -2  nT  occurred.  For  example,  for  a 
given  hourly  measurement  of  B ,  positive  (either  weak  or 
strong),  the  probability  of  encountering  a  B.  of  less  than  -2  nT 
within  the  past  6  hours  is  around  63%.  It  is  evident  from  the 
plot  that  there  is  no  significant  difference  hi  B.  history  for  the 
two  B.  positive  bins  used  here.  Thus  the  increase  in  hemi¬ 
spheric  fluxes  for  strong  B.  positive  cannot  be  attributed  to  a 
residual  effect  due  to  previous  B.  negative  periods. 

6.2.  Correlation  of  Hemispheric  Fluxes 
and  the  Solar  Wind  Pressure 

The  values  of  the  solar  wind  dynamic  pressure  {mpV;w) 
listed  in  Table  2  follow  a  similar  trend  with  respect  to  B.,  as 


Bz  NORTH  PERSISTENCE  (hours) 

Fig.  9.  A  plot  of  the  relative  proximity  of  preceding  it-  soulh 
intervals  to  both  moderate  and  strong  B:  north  intervals. 


do  the  ion  hemispheric  fluxes.  That  is,  the  pressure  ap¬ 
proaches  a  minimum  for  small  to  moderately  large  B, 
positive.  One  difference  is  that  the  pressure  is  more  sym¬ 
metric  about  its  minimum,  with  more  nearly  equal  pressure 
for  both  strong  B.  south  and  strong  B,  north.  The  increase 
in  dynamic  pressure,  and  subsequent  compression  of  the 
magnetosphere,  could  partly  explain  the  increase  in  total 
precipitating  number  and  energy  particle  flux.  This  possibil¬ 
ity  is  supported  by  the  work  of  Scurry  and  Russell  [1989]. 
They  analyzed  the  activity  index  Am  and  the  solar  wind 
pressure  during  periods  of  northward  B.  and  concluded  that 
the  apparent  correlation  between  geomagnetic  activity  and 
increasingly  northward  B.  is  largely  a  side  effect  of  increas¬ 
ing  dynamic  pressure  during  periods  of  large  B .. 

6.3.  Fits  to  Data 

We  have  shown  that  for  both  electrons  and  ions,  the 
variation  of  both  the  average  hemispheric  particle  energy 
flux  and  the  particle  number  flux  with  the  solar  wind  velocity 
V’w  and  IMF  B.  component  can  be  represented  by  equation 
(3)  or  (4).  A  large  body  of  work  exists  on  the  nature  of  the 
solar  wind-magnetosphere  interactions.  Many  of  these  have 
been  investigations  of  the  relationship  between  some  index 
reflecting  the  approximate  level  of  geomagnetic  activity  and 
some  combination  of  parameters  in  the  solar  wind.  Typi¬ 
cally,  indices  such  as  AE,  AL,  Kp,  and  Ap  have  been 
investigated  as  a  function  of  ViWBz,  V%WBS,  V;WBS,  and  e 
[Baker  et  al.,  1983,  and  references  therein:  Bargatze  et  al., 
1985].  In  such  studies,  B;  is  the  average  value  of  the 
north-south  component  of  the  IMF,  while  Bs  is  the  value 
averaged  only  when  B.  is  southward.  The  parameter  e  = 
ViVIB2  sin4  (6/2)/|j ,  where  0  is  given  by  tan-1  (ByIBz)  or  by 
this  angle’s  supplement  for  B.  <  0  and  lQ  is  a  constant.  In 
many  of  these  studies,  the  assumption  has  been  that  the 
magnetosphere  operates  as  a  half-wave  rectifier  with  activity 
increasing  steadily  for  increasingly  negative  B.,  but  exhibit¬ 
ing  negligible  activity  for  increasingly  positive  B.  [Burton  et 
al.,  1975]. 

The  studies  cited  above  have  shown  in  general  that 
VsvtB.,  VSWB,,  V2wBs,  and  e  all  provide  approximately  the 
same  level  of  correlation  with  the  various  geomagnetic 
activity  indices.  Our  function  is  somewhat  similar  to  the 
epsilon  function  in  that  it  has  a  linear  dependence  on  Vsw 
and  a  quadratic  dependence  on  Bz.  Unlike  the  epsilon 
parameter,  however,  our  functional  form  shows  no  i>.:ong 
gating  with  the  IMF  polarity,  as  would  be  required  for  a 
half-wave  rectifier  response.  Rather,  we  show  a  minimum  in 
particle  precipitation  for  B.  weakly  to  moderately  positive, 
with  an  increase  for  strong  B.  positive.  This  is  consistent 
with  studies  of  other  low-altitude  signatures  used  to  monitor 
magnetospheric  activity.  For  B.  increasing  from  weak  to 
strong  positive  values,  Reiff  and  Luhmann  [1986]  have 
pointed  to  studies  suggesting  an  increase  in  the  cross-polar 
cap  potential,  Gussenhoven  [1982]  has  found  an  increase  in 
polar  cap  arc  activity,  Zanetti  and  Poteinra  [1986]  point  to 
the  establishment  and  enhancement  of  “NBZ”  currents  in 
the  polar  region,  and  Makita  et  al.  [1988]  have  shown  an 
enhancement  of  the  soft  electron  boundary  plasma  at  very 
high  latitudes.  These  studies  all  show  that  the  magneto¬ 
sphere  becomes  energized  through  mechanisms  driven  by  B. 
positive  as  well  as  by  Bz  negative. 
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TABLE  5.  Efficiency  of  V5W  for  Bz  South  and  North 


Vsw,  km/s 

B  £flvx 

,  GW 

"nflux*  <s  sr)" 

l 

Electron 

Ion 

Electron, 

o 

X 

Ion,  xlO23 

Bzs 

Bzn 

Bzn 

Bzs 

BZn 

Bzs 

Bzn 

346 

18.6 

6.5 

2.0 

1.2 

2.5 

1.8 

4.6 

4.0 

677 

50.7 

29.1 

3.0 

2.4 

4.3 

3.3 

5.1 

4.2 

Ratios 

2.7 

4.5 

1.5 

2.0 

1.7 

1.8 

1.1 

1.1 

6.4.  Efficiency  of  Vsw  in  Controlling 
Energy  Transfer 

Rossberg  [1984]  showed  that  during  periods  of  Bz  north, 
AE  increases  linearly  with  Vsw.  In  his  subsequent  work 
[Rossberg,  1989],  he  examined  the  variation  of  the  magnetic 
activity  index  AL  for  both  Bz  south  and  Bz  north  during 
intervals  of  both  low  Vsw  (<400  km/s,  with  an  average  of  358 
km/s)  and  high  Vsw  (>700  km/s,  with  an  average  of  740 
km/s).  For  periods  of  both  8,  south  (BzS)  and  Bz  north 
( BzN ),  he  determined  the  ratios  of  AL  for  high  and  low  Vsw. 
From  these  ratios,  he  concluded  that  the  solar  wind  speed 
dependent  processes  involved  in  solar  wind-magnetosphere 
coupling  are  twice  as  efficient  in  their  transfer  of  energy 
during  periods  of  Bz  north  than  periods  of  Bz  south. 

We  assume  that  the  hemispheric  energy  flux  of  precipitat¬ 
ing  particles  is  at  least  as  good  as  the  AL  index  in  providing 
a  measure  of  the  efficiency  with  which  Vsw  contributes  to  the 
energy  transfer  from  the  solar  wind  to  the  magnetosphere. 
Using  this  assumption,  we  follow  an  analysis  similar  to 
Rossberg’s,  substituting  the  hemispheric  energy  flux  for  AL. 
To  provide  the  best  comparison  to  Rossberg’s  results,  we 
combined  our  data  to  determine  the  average  value  of  the 
electron  and  ion  HEhux  for  BzS  and  BzN  for  the  lowest  (346 
km/s)  and  highest  (677  km/s)  Vsw  cases.  These  data,  as  well 
as  the  corresponding  fINnm  data,  are  listed  in  Table  5.  For 
all  columns  the  flux  ratios  for  the  high-velocity  case  to  the 
low-velocity  case  are  given  in  the  last  row.  Rossberg  found 
that  for  BzS  (BzN),  a  doubling  of  Vsw  leads  to  an  increase  by 


a  factor  of  2.8  (5.9)  in  the  average  AL  index.  Our  study  gives 
a  similar  result  in  the  average  hemispheric  electron  flux  for 
BzS  (BzN).  A  doubling  of  Vsw  leads  to  an  increase  by  a  factor 
of  2.7  (4.5).  Thus  both  studies  show  that  an  increase  in  Vsw 
is  most  efficient  in  increasing  magnetospheric  activity  during 
northward  Bz .  As  is  indicated  in  Table  5,  the  ion  hemispheric 
energy  flux  exhibits  a  similar  if  somewhat  less  pronounced 
variation.  Although  there  is  a  clear  asymmetry  in  the  effi¬ 
ciency  of  Vsw  in  transferring  energy  to  the  magnetosphere 
under  the  two  different  Bz  polarities,  such  is  not  the  case  for 
the  transfer  of  mass  (number  flux).  There  are  no  significant 
differences  for  Bz  north  and  south  in  corresponding  ratios  for 
the  average  hemispheric  number  flux. 

6.5.  Comparisons  With  Other  Studies 

The  majority  of  the  statistical  studies  mapping  out  the 
region  of  precipitating  particles  have  been  ordered  by  the 
magnetic  activity  index  Kp.  For  electrons,  these  include 
Wallis  and  Budzinski  [1981],  Spiro  et  al.  [1982],  Hardy  et  al. 
[1985],  and  Foster  et  al.  [1986];  and  for  ions,  Hardy  et  al. 
[1989].  Brautigam  et  al.  [1988]  compared  the  electron  hemi¬ 
spheric  energy  flux  from  these  studies  and  found  excellent 
agreement. 

As  an  indirect  comparison  of  the  hemispheric  fluxes 
derived  from  our  Kp  statistical  studies  (papers  1  and  2)  and 
those  from  the  study  presented  here,  we  plot  electron  and 
ion  hemispheric  flux  versus  Kp,  using  the  Bz:ViVI  bin- 


KP  KP 


Fig.  10.  Hemispheric  number  flux  (Hnhvx)  from  both  the  /wp-ordereci  study  and  the  current  sludy,  plotted  as  a 

function  of  Kp.  (Left)  Electrons.  (Right)  Ions. 
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Fig.  11.  Hemispheric  energy  flux  (//£pux)  from  both  the  ATp-ordered  study  and  the  current  study,  plotted  as  a 

function  of  Kp.  (Left)  Electrons.  (Right)  Ions. 


averaged  Kp  values  listed  in  Table  2.  These  are  shown  in 
Figure  10  (#/vHUX  versus  Kp)  and  Figure  11  (//£rux  versus 
Kp).  In  these  plots  the  seven  flux  values  from  the  earlier  Kp 
studies  have  been  plotted  as  a  solid  line,  and  the  26  flux 
values  from  the  current  study  are  marked  by  a  symbol 
according  to  their  B.  bin  (independent  of  their  Vsw  bin).  The 
results  from  the  Kp  statistical  maps  are  tabulated  in  Table  6. 
Generally,  there  is  good  agreement  between  the  two  sorts. 
The  largest  discrepancies  occur  for  the  higher  negative  Bz 
cases  where  the  values  consistently  tend  to  lie  above  the 
values  for  the  previous  Kp  sort.  This  is  most  marked  for  the 
electron  HEum  case.  This  difference  reflects  the  mixing  of 
B .  polarities  in  the  Kp  sort  as  we  showed  for  Kp  =  4  and  6 
in  Figure  8.  For  the  higher  Kp  level  maps  there  is  significant 
mixing  of  Bz  positive  with  the  more  dominant  Bz  negative. 
Since  for  Bz  positive  the  general  level  of  precipitation  is 
lower,  to  include  these  Bz  positive  periods  in  the  Kp  study 
decreases  the  average  hemispheric  flux  as  compared  to  the 
current  study.  The  disparity  in  results  is  greatest  for  the 
electron  hemispheric  energy  flux  because  the  electron 
//£f|(lx  shows  the  greatest  response  to  Bz  negative  and  thus 
most  greatly  amplifies  the  problem  of  mixing  data  under 
different  Bz  polarities. 

Foster  et  al.  [1986]  also  reported  the  electron  hemispheric 
energy  flux  for  two  Bz  intervals:  >+3  nT  and  <-3  nT.  They 
analyzed  NOAA/TIROS  satellite  data  and  found  an  increase 
in  hemispheric  electron  energy  flux  from  14  GW  for  their  Bz 


positive  interval  to  32  GW  for  their  Bz  negative  interval.  For 
a  direct  comparison,  we  merged  all  our  Vsw  bins  together 
and  computed  electron  //£|iux  (from  unextrapolated  average 
spectra)  for  the  two  bins  with  Bz  intervals  of  + 10  to  +3  nT 
and  -3  to  - 10  nT.  We  found  the  energy  fluxes  for  these  Bz 
positive  and  B.  negative  intervals  were  11  and  53  GW, 
respectively.  Foster  et  al.’s  wattage  forfl,  positive  is  slightly 
higher  (xl.3)  than  ours,  whereas  their  wattage  for  Bz 
negative  is  more  significantly  lower  (x0.6).  One  source  of 
discrepancy  which  could  explain  their  lower  energy  input  for 
Bz  <  -3  nT  is  that  the  NOAA/TIROS  detectors  have  a 
smaller  enerjt\  range,  with  an  upper  threshold  of  20  keV 
versus  our  detectors'  threshold  of  30  keV. 

7.  Summary  and  Conclusions 

Using  a  large  data  set  of  precipitating  electron  and  ion 
measurements,  we  have  been  able  to  determine  the  variation 
in  the  hemispheric  power  input  and  number  flux  due  to  these 
particles  as  a  function  of  polarity  and  magnitude  of  the  IMF 
B .  component  and  as  a  function  of  the  magnitude  of  the  solar 
wind  velocity.  Both  the  hemispheric  energy  flux  and  hemi¬ 
spheric  number  flux  are  found  to  vary  in  a  simple  and 
consistent  manner  with  both  B.  and  Vsw.  For  both  electron 
and  ion  fluxes  the  variation  with  B.  (for  a  fixed  Vsw)  tends  to 
a  minimum  for  weak  to  moderately  strong  Bz  positive  and 
increases  at  a  rate  faster  than  linear  for  values  both  above 


TABLE  6.  Hemispheric  Fluxes  for  Kp  Study  (Electron  and  Ion) 


Map 

Number 

Kp 

Electron 
«e„ux-  GW 

Electron 
xlO25  (s  sr)_l 

Ion 

"W-GW 

Ion 

tf/W  xlO2’’  (s  sr)"1 

1 

0 

4.68 

1.53 

0.73  (0.67) 

3.62  (3.46) 

2 

1 

8.14 

1.85 

1.38  (1.26) 

3.80  (3.76) 

3 

2 

14.90 

2.35 

2.03  (1.61) 

4.54  (4.25) 

4 

3 

25.90 

3.08 

2.92  (2.06) 

5.53  (5.05) 

5 

4 

33.60 

4.05 

4.10(2.43) 

7.10(6.00) 

6 

5 

43.20 

5  83 

5.47  (2.75) 

8.74  (7.10) 

7 

26 

/7.60 

10.20 

10.00  (5.69) 

15.00  (12.90) 

Unextrapolated  values  are  listed  in  parentheses. 
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and  below  the  minimum.  This  effect  is  most  marked  for  the 
ions.  The  variation  in  the  fluxes  with  Vsw  tends  to  be  linear 
with  the  steepest  slopes  for  the  electron  hemispheric  energy 
flux.  The  variation  with  Bz  and  Vsw  can  be  fit  well  for  all  four 
quantities  by  an  equation  either  of  the  form 

ABV  Vsw)=a(Bz-b)2  +  cVsw  +  d 

or 

ABZ,  Vsw)  =  [a(Bz  -  b)2  +  l][cVsw  +  d]. 

Either  equation  indicates  that  no  variation  in  Bz  higher  than 
quadratic  and  no  variation  in  Vsw  higher  than  linear  is 
required.  In  fact,  the  first  form  shows  B?  and  Vsw  to  be 
independent.  This  is  significantly  different  from  either  the 
half-wave  rectifier  or  epsilon  function  response. 

In  addition,  we  find  that  the  ions  make  the  largest  contri¬ 
bution  relative  to  the  electrons  for  the  lowest  velocity  and 
the  highest  value  of  Bz  positive.  The  relative  contribution 
decreases  sharply  for  decreasing  values  of  Bz  positive  and 
for  increasing  velocity.  For  Bz  negative  for  a  fixed  velocity 
the  relative  contribution  is  approximately  constant.  The 
relative  contribution  of  the  ion  to  the  electron  number  flux 
shows  a  similar  variation.  Abo  significant  is  that  Vsw  is  most 
efficient  in  increasing  energy  flux  during  periods  of  Bz 
positive,  and  more  so  for  electrons  than  ions. 

The  electron  and  ion  total  hemispheric  number  and  energy 
flux  represent  only  one  product  from  the  large  data  bases 
which  have  been  compiled.  A  subsequent  study  utilizing  the 
same  data  bases  will  examine  the  dependence  of  specific 
auroral  regions  on  the  IMF  Bz  and  Vsw  in  terms  of  individual 
spectra  and  auroral  boundaries.  This  more  detailed  analysis 
is  required  before  one  can  comment  further  on  several  of  the 
issues  raised  in  this  paper. 


Appendix:  Fitting  Procedure 

The  similarity  in  the  variation  of  the  four  quantities  with 
Bz  and  Vsw  suggests  that  the  variation  can  be  fit  by  a 
relatively  simple  functional  form.  Two  different  approaches 
were  taken  to  this  fitting.  In  the  first  approach,  we  found  that 
the  variation  with  Bz  in  all  four  quantities,  for  a  fixed  Vsw, 
could  be  well  fit  by  the  general  parabolic  equation 

f{Bz)  =  axB2z  +a2Bz  +  ay 

Using  this  fact,  we  fit  separately  the  Bz  variation  of  each  of 
the  four  hemispherically  integrated  quantities  for  each  dif¬ 
ferent  velocity  to  an  equation  of  this  form.  The  resulting  fit 
parameters  a  |  and  a2,  for  each  of  the  four  quantities,  were 
found  to  have  no  significant  V5W  dependence  and  were  thus 
averaged  over  the  five  velocity  cases.  Parameter  a3  was 
found  to  vary  linearly  with  Vsw  and  could  be  fit  to  the  linear 
equation 

ai  =  a3^sw  +  fl4- 

Combining  those  two  relations  gives  an  equation  of  the  form 
f(B j,  V5W)  —  Q\Bz  +  Oj Bz  "b  ®3^sw  "b 
This  equation  may  be  expressed  as 


J{BZ,  Vsw)  =  a(Bz  -  b)2  +  cVsw  +  d.  (Al) 

This  relation  shows  that  the  data  can  be  simply  described  as 
a  family  of  parabolas  whose  minima  are  located  at  Bz  =  b 
and  are  offset  by  a  linear  function  of  Vsw. 

The  second  method  is  also  based  on  the  observation  that 
the  fluxes  for  each  hemispheric  quantity  may  be  fit  to  a 
parabola.  In  this  method,  the  fluxes  as  a  function  of  Bz  were 
normalized  to  the  minimum  flux  value  for  a  fixed  Vsw.  Since 
the  normalized  values  varied  only  slightly  with  Vsw  for  a 
fixed  Bz ,  the  normalized  values  were  then  averaged  over  the 
different  solar  wind  velocities  and  the  average  fit  to  the 
equation 

f{Bz)lf{Vsw)min=  a](Bz  -  a2)2  +  a}. 

As  in  the  first  procedure,  the  variation  in/fVsw)mjn  can  be 
approximated  by  the  linear  equation 

■/0/$w)min  —  a4^sw  "b  <*5. 

Combining  the  above  gives  the  equation 

flBz,  Vsw)  =  [a{{Bz  -  a2)2  +  a3][a4VSw  +  05]- 

Factoring  the  a3  out  of  the  first  bracketed  expression  and 
into  the  second,  the  equation  may  be  simplified  to 

ABZ,  Vsw)  =  [a(Bz  -  b)2  +  l][cVsw  +  dl  (A2) 

Although  the  final  functional  forms  for  both  methods, 
given  as  equations  (Al)  and  (A2),  have  four  free  parameters 
each,  the  functional  forms  differ.  Equation  (Al)  is  the 
simpler,  containing  only  the  terms  B?,  Bz,  and  Vsw.  Equa¬ 
tion  (A2)  contains  the  additional  cross  terms,  Viv/Bz  and 
Viv/B2.  Because  there  are  four  free  parameters  in  each 
method,  the  quality  of  fit  should  be  approximately  the  same. 
Nevertheless,  we  can  attempt  to  judge  which  method  is 
better  by  investigating  the  goodness  of  fit  parametrically  and 
point  by  point. 

The  four  constants,  a  through  d,  for  the  two  methods  of 
fitting  are  listed  in  Table  4.  Two  constants  have  the  same 
meaning  in  each  method  and  can  be  directly  compared.  The 
constant  b  equals  the  value  of  Bz  that  gives  the  minimum 
flux  for  a  given  V5W .  The  constant  d  is  the  residual  flux  when 
the  effects  of  B,  and  Vsw  are  set  to  zero.  Or  alternately 
interpreted,  -die  is  the  solar  wind  velocity  that  is  required 
before  there  will  be  a  positive  flux  into  the  hemisphere  when 
Bz  =  b. 

Comparing  the  values  of  b  for  the  different  types  of  fluxes, 
we  find  excellent  overall  agreement  between  the  two  meth¬ 
ods.  The  worst  comparison  is  for  the  case  of  electron  energy 
flux  where  both  methods  show  that  the  minimum  flux  is 
obtained  for  large  positive  B.  and  the  fit  is  done  with  only 
one  half  of  the  parabola.  It  is,  however,  interesting  to  note 
that  while  there  is  agreement  between  the  two  methods,  the 
value  of  B.  for  minimum  flux  varies  greatly  for  the  type  of 
flux  considered.  It  is  largest  for  the  electron  with  an 

average  of  3.9  ±  0.6  nT  between  the  two  methods;  has  an 
average  value  of  2.1  ±  0.2  nT  for  the  electron  //W|LUX  and  ion 
and  is  smallest  for  the  ion  Hnhvx>  where  both 
methods  give  ~0.95  nT.  The  agreement  in  both  value  and 
trend  for  the  Bz  that  minimizes  each  flux  does  not  help  us  to 
choose  between  the  methods. 
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Comparing  the  constants  d  for  the  two  methods  indicates 
a  problem  whenever  the  values  are  negative,  since  for 
negative  values  some  minimum  V sw  must  be  reached  before 
the  flux  will  be  positive  when  B,  is  near  b.  Negative  d  values 
occur  for  the  energy  fluxes  of  both  species.  The  minimum 
Vsw  value  is  equal  to  -die.  For  the  electron  energy  flux  this 
value  is  372  km/s  and  211  km/s  for  methods  1  and  2, 
respectively.  The  former  value  is  much  too  high  and  indi¬ 
cates  that  for  low  solar  wind  velocities  method  1  will  not  give 
good  values.  This  is  the  case  when  the  values  calculated 
from  the  functional  form  in  equation  (Al)  are  compared  with 
data  points  for  the  lowest  velocity  bin.  For  method  2  the 
minimum  Vsvi  for  the  electron  energy  flux  is  lower  than 
measured  velocities  of  the  solar  wind,  and  comparison  of 
calculated  and  measured  values  shows  that  there  is  no 
abnormal  behavior  in  the  low-velocity  bin.  This  same  argu¬ 
ment  holds  for  the  minimum  velocities  for  the  ion  energy  flux 
which  are  95  and  80  km/s  for  methods  1  and  2,  respectively. 

Acknowledgments.  The  Editor  thanks  R.  D.  Zwickl  and  another 
referee  for  their  assistance  in  evaluating  this  paper. 
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